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It is expected that black holes are formed dynamically under gravitational collapses and approach 
stationary states. In this paper, we show that the asymptotic Killing vector at late time should 
exist on the horizon and then that it can be extended outside black holes under the assumption of 
the analyticity of spacetimes. This fact implies that if there is another asymptotic Killing vector 
which becomes a stationary Killing at a far region and spacelike in the "ergoregion," the rotating 
black holes may have the asymptotically axisymmetric Killing vector at late time. Thus, we may 
expect that the asymptotic rigidity of the black holes holds. 
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I. INTRODUCTION 

Black holes in our Universe are expected to be formed under gravitational collapses, and to finally approach 
stationary and vacuum states by radiating and absorbing energy, momentum, and angular momentum. Then the 
uniqueness theorem [l[ guarantees that the black hole candidates in our Universe are the Kerr black holes. The key 
ingredient for the proof of the uniqueness theorem is the rigidity theorem The rigidity theorem shows that 

the stationary rotating black holes have axisymmetric Killing vectors. The outline of the proof is as follows: The 
^ ] stationarity of spacetime implies that there are no gravitational waves around a black hole. Then the expansion and 
shear of the event horizon must vanish by virtue of the Raychaudhuri equation and stationarity. Using the Einstein 
5^ I equations, then, we can find that the null geodesic generator of the event horizon is a Killing vector. If the black 
bJ} hole is rotating, this new Killing vector may deviate from the stationary Killing vector which becomes spacelike in 
the ergoregion. This means that the stationary black holes might have two Killing vectors, that is, the stationary and 
axisymmetric Killing vectors. Hence the stationary black holes should rotate rigidly. 

However, the latc-timc phase of black holes produced by gravitational collapses would not be exactly stationary 
but nearly stationary. "Nearly stationary" means that the black holes arc surrounded by the gravitational waves at 
late time. Then we cannot apply the rigidity theorem to such black holes because the expansion and shear of the null 
geodesic generator on the event horizon do not vanish due to the presence of the gravitational waves on the event 
horizon. Note that the late-time behaviors of the perturbations around the Schwarzschild and Kerr black holes were 
examined and it was shown that the perturbations decay at late time both on the horizon and null infinity at the 
same rate (for examples, see Refs. 043)- the dynamical processes in gravitational collapses, however, it is quite 
nontrivial whether the formed black hole approaches the Kerr black hole. In this paper, we show that the asymptotic 
Killing vector, which will asymptotically approach a Killing vector at late time, should exist on the horizon without 
assuming any symmetries and then that it can be extended outside of the event horizon using the Einstein equations. 
This indicates that if there is another asymptotic Killing vector which will be an asymptotically stationary Killing 
vector at a far region, the rotating black hole may have the asymptotically axisymmetric Killing vector at late time. 

This paper is organized as follows. In the next section we provide the Einstein equations near the future event 
horizons. In Sec. HI, we investigate the late-time behaviors of the metric on the event horizon and find the late-time 
symmetry. Then, under the assumption of the analyticity, we show that this late-time symmetry can be extended 
outside of the event horizon using the Einstein equation. In Sec. IV, we provide a summary and discussion. In 
Appendix A we perform the (n -f- 1) decompositions of the Einstein equation. Using them, we provide the explicit 
form of the Einstein equation in the current form of the metric in Appendix B. Almost a similar formulation is 
developed in the study on null infinity In Appendix C, we present the details of the discussion associated with 
gauge freedom. 



II. BONDI-LIKE COORDINATE AND EINSTEIN EQUATIONS 



In this section, introducing the Bondi-likc coordinate near the event horizon, we investigate the initial value problem. 
For the details of these derivations, see Appendixes A and B. 
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A. Bondi-like coordinate near event horizon 



We consider a dynamical black hole and investigate the late-time behavior of the near horizon geometry at late 
time. The black hole is defined as the region which is not contained in J~(j^"*") where J^"*" is the future null infinity. 
J~{S) is the chronological past connecting a set S by causal curves from S. The event horizon is the boundary of 
the black hole and it is a null hypcrsurface. See Ref. 0] for the precise definitions. Then, we introduce the Bondi-like 
(Gaussian null) coordinates = {u, r, a;°) near the event horizon as 

ds^ ^ cjABdx^dx^ = -Adu^ + 2dudr + habidx" -f U°du){dx'' + U^du), (1) 

where w is a time coordinate. In this coordinate, the horizon position is taken to be r = 0. We assume that a cross 
section of the event horizon is compact in the u = constant hypcrsurface and its topology is S^. arc coordinates 
on S'^ . Since the event horizon is a null hypcrsurface, g^u must vanish on the event horizon and I = d/du is the null 
geodesic generator on the event horizon. Furthermore we can choose the coordinate u so that = on the event 
horizon. Then, we have 

A = 0, (2) 

and 

[/° = 0, (3) 

where = means the evaluation on the event horizon r = 0. 

To solve the vacuum Einstein equations, we formulate the initial value problem in the Bondi coordinates. For the 
convenience of our study in the following sections, wc solve the Einstein equations as the evolution equations in the 
direction of r. Thus, the initial value of the metric should be set on the r = constant surface. In this paper wc take 
the event horizon r = as the initial surface. 

The evolution equations are given by 

Rrr = -\{\0gh)" --^h'^'^h'^hab)' [Kd)' = 0, (4) 

RrBh'^'' = ^f/"" + \h''''KcU^' + \{\og hyu'^' + h^'D'ih',^ - /ifc,(log h)'] = (5) 



and 



RABh\h^fj 



-\AKh - \Ahab - {hab)' + ^''^Rab + ^h^'^Kcha + ^h^'' (Kckd + K^c) + ^uKb 

- ^KciDbU' + D'U,) - ^h,^{DaU- + D'Ua) - ^hachMU^'u"' - i[(log/i) - 2DaU'']h'^, 

- ^{\0g h)'iAh'^, + Kb - DaUb - DhUa) + \{hbcDaU'' + hacDbU^') = 0, (6) 

where the prime and dot denote the r and u derivative, respectively. Da is a covariant derivative with hab and 
h ~ det hab- Also, ^'^^Rab is the Ricci tensor with respect to hab- The evolutions of the metric functions A, t/° and 
hab are determined by Eqs. ([S]) and ([H]) completely. Note that the trace part of Eq. ^ gives us 

A'(log/i)' = 2'^''^ R ~ ^[{\og h)']^ -[{log h)-2DaU'']{\ogh)' ^habU'^'u'' 

- A{\og h)" - 2(log h)' + 2{DaU'')' + U^Da{\og h)', (7) 

which is used as the evolution equation for A. 

The other components of the Einstein equations, R^u = 0, Rur = and Rua ~ 0, are related to the above evolution 
equations by the Bianchi identities. Therefore, once they are satisfied at the initial surface r = 0, we do not need to 
solve them any more. In fact, if the evolution equations are satisfied, the Bianchi identities lead to 

{VhRuaY = VhR ur 1 

{VhRraY = -VhDaR\ - {Vh)Rur, 
{loghYRur = 0, 
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where R^u = Ruu + ARur — U'^Rua and = —V^Rur + h°^Rab- Thus, the evolution equations guarantee that 
Ruu = 0, Rur = and R^a = will always be satisfied at r 7^ if = 0, Rur = and Rua = at the initial surface 
r = 0. For convenience, at r = we will use the following constraint equations: 

Ruu= - \i\0gh) + y (log/l) - -^h-'^h^'-habhcd = 0, (9) 

and 

- ^(f/")' - |(log /i)f/"' - ^h-^U^'kbc - ^&'A' + ^h-'D^hc - i^"(log h) = (10) 

[see Eqs. (|A30[) and (|A31|) ]. Moreover, since A should vanish on the initial surface r = [see Eq. ([5])], the evolution 
equations Eq. ([6]) become singular on r = 0. Analyticity of hab on r = gives us the regularity condition 

- iKchbdU^'u"' -\{\ogh)K, 

- ^{\ogh)' Kb + \{hbcDaU''' + KcDbU''') = 0. (11) 

If we give hab\r=Q on r = 0, we can determine /i^j|r=o, C^"'|r=o, and A'|r=o by solving Eqs. ([TUl) and PT|) . As 
a result, we will solve the evolution equations (|4]), (O and ([6]) using the above initial values on r = 0. 



B. Some explicit solutions near event horizon 

In this subsection, it is shown that we explicitly solve the constraint equations and evolution equations near the 
event horizon using power series expansion around r = 0. To do this, we expand the metric functions with r near the 
event horizon as 

A = rA(i) +^r*A(*\ (12) 

i>2 



Ija ^ rC/^i)" + ^r't/W", (13) 



i>2 

and 



ha, = h^+rl^+Y.^h^^ (14) 



i>2 



where from the gauge conditions Eqs. ^ and ([3]), the expansions of A and U"" start from the first order of r. In the 
following, the tensor index of h^^l and t/^*-''^ is raised and lowered by h!fb- trace and traceless parts are also taken 

First let us solve the constraint equations and regularity conditions in order to determine initial values. On the 
initial surface r = 0, h^^^ , A^^^ and U^^^"" should be set on r = as initial values. The constraint equations for 
these initial values are Ruu = and R^"' = 0. Now i?„„ = [Eq. ^] is rewritten as 

/i(o) „ + = 0, (15) 

where = deth^^^. Thus, A^^'^ should be given to satisfy this equation for given h^^^ . Also, i?'"° = [Eq. pU)) ] 
becomes 

^ + h'-°^'"'D''h[°J - £i'^(log/i(o)) - /i(°)'"=C/(i)''/i[,°^ - ic/(i)'^(log/i(o)), (16) 
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where Da is the covariant derivative with respect to h^^b ■ Then the initial value f/^^)" is given to satisfy the above. 
The regularity condition Eq. PT|) becomes 

- \iloghio))h^^ - Ih^'^^'hi'Jh^J + liDaUi'^ + D,ul^y) = 0, (17) 

where ^'^^"^^Rab is the Ricci tensor of h!fb- obtain /i^^' satisfying this equation. Hence, if we give h^^^ on the initial 

surface, the constraint equations and the regularity condition yield h^^b ^ -^^^^ ^'^'^ U^^^"'. 
Next we solve the evolution equations. Equation (j4]) becomes near the event horizon as 

Rrr = -h^'^^'h^^^ + \{h^abf + 0{r) = 0. (18) 

This equation gives us the trace part of as 

h^'^-'h^a^ = \ih'^f- (19) 

The evolution equation of [/" as Rrsh"^^ ~ [Eq. ([5])] can be expanded as 
Then is given by 

Ui2)a ^ _ l;,(0)ac/^(l)f^(l)fc _ 1 c;(l)a/^(0)bc^(l) _ }_J^il)ac j^b (^^W _ /^^ /.W^'^^l) ) . (21) 

In the same way, expanding the evolution equations Eq. ([6]) near the event horizon, we can obtain the traceless part 
h^^ab) ^'^'^ ^^"^^ ™ terms of h^ab ^^ab ^^'^ A^^\ Note that A^'^^ is given by the trace part of Eq. ([6]), namely 

Eq. ([7]). However, we do not provide its explicit form because its form is very cumbersome. 

Hence, we can determine h^^b ^ C/*-^'", and A^^' using the evolution equations. To determine the higher order 
quantities, h'"^^, U^^^"', and A^*' (z > 2), wc have to repeat the same procedure. 



III. LATE-TIME SYMMETRY ON AND NEAR EVENT HORIZON 



In this section, we show that there is late-time symmetry on the event horizon. Then we will extend it outside of 
black hole regions. 



A. Late-time behaviors on event horizon 



To investigate late-time behaviors of the event horizon, we introduce the expansion and shear of the null geodesic 
generator I = d/du of the event horizon. The expansion 9 and shear aab are defined as 

where Uab is the traceless part of h^^^ with respect to h^^^ . Then we can rewrite Eq. one of the constraint 
equations, using and Uab as 

e-]^A^^'^e = -aab(7''\ (23) 
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We can regard as a surface gravity of the black hole with respect to the time coordinate u because 



4(1) 

l^V^lB^^lB (24) 



cMT^^u'l (25) 



holds. Using the affine parameter w defined by 

d 
du 

we can obtain the Raychaudhuri equation 

dw6(^ui) = -^Sfw) - <^{w)ab<^{iu), (26) 

where 9^^j and (T(^^jab expansion and shear with respect to w. We can see the relation between 6, aab and ^(u,), 

0'(to)a& as 

9 = e'f^jexpM -^du'j , aab = cr(w)abexpn -^du'j. (27) 



Here we remember that the area law of the event horizon holds for spacetimes satisfying the null energy condition, 
that is, 6* > and ^(u,) > 0. Since (Ji^^-^ab^l^) ^ 0, Eq. implies the inequality 

+ \el^^ < 0. (28) 

Then the integration over w gives us 

V^'(o) + {w- wo)/2 

where we used the fact of 0(o) = 6'(^)(w = Wq) > 0. In addition, Eq. (|26p shows that the shear tT(„)ah should also 
vanish as w — > oo. This is shown as a part of the proof of another theorem [lo| . From now on, we assume that 
w — !■ oo corresponds to it — oo. Then, we see that 9 and (Tab should also vanish as w — >■ oo. It is natural to assume 
that the cross section of the event horizon is compact. Then the vanishing of the expansion implies that the horizon 
area approaches a constant and finite value. 

Altogether we see the behavior of the metric at late time as 

'Ci.gAs I horizon (u^oo). (30) 

Here we impose the following decaying condition on the event horizon for the metric: 

A5AB = 0(^-1), (31) 

which explicitly means h!"^^ ~ 0{u~"). This equation means that the null geodesic generator of the event horizon 
I should be an asymptotic Killing vector at late time (u — > oo). Thus there is a late-time symmetry on the event 
horizon. 

Since we consider the dynamics only near the event horizon, we cannot determine the decaying rate of the metric. 
However, the details of the decaying properties are not important for our argument, that is, our result does not depend 
on n. Note that it will be determined by the boundary conditions at a far region from the black holes as in Refs. 



It should be remembered that the horizon which satisfies this condition is called a slowly evolving horizon in Refs. 
[Ill, [13] ■ If the right-hand side of Eq. (|3ip vanishes, the event horizon will be identical to the isolated horizon [isj . 

B. Extension of late-time symmetry 

The purpose of this subsection is to show that the decaying condition of Eq. ([5T|) can be extended outside of the 
event horizon as 

CigAB =o(^). (32) 
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In the following we assume the analyticity of gAB- 

Under the presence of the analyticity of spacetimes, the above is equivalent with 



{Cr^y^CigAB^oQ^y (33) 

where n = d/dr and jti = 0, 1, 2, • • • . Note that "=" means the evaluation on the event horizon (r = 0) again. 
First we will show the m = 1 case: 

CnClQAB = O . (34) 

Substituting the explicit form of gAB to the above, we rewrite Eq. p4p as 

- CnCl{A-UaU''){du)A{du)B+2CnClUa{du)(A{dx'')B)+CMab{dx'')A{dx^^ ■ (35) 

Using Eqs. (fT2|) . (|T3l) and ([T4)) . the above will be equivalent with 

=o(^\, (36) 



u" 



(37) 



and 



^^h^a^ = O ( ) ■ (38) 
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Let us examine these conditions. First, we focus on Eq. ([36]). As a result, using the gauge freedom, we can show that 
the slightly strong condition, = O (l/u"+^), holds. To see this, we decompose into the M-independent 

term and others as 

= a'^^\x'')+A'^^\u,x''). (39) 

As shown in Ref. we can choose the gauge so that is a constant. Furthermore, using the residual gauge, 

we can take i'l' = 0(l/ii"). Then, using the gauge freedom in our coordinate, we can make A^^'^ satisfy stronger 
decaying property as 



^^A^'^ - O {-^] ■ (40) 



For the details, see Appendix C. Of course, Eq. ([55)) is satisfied. 
Now 11^^'^°' should satisfy the constraint equation of Eq. ([16p as 



= -D-^i^i) + h^°'>'"'D''hfJ - D°(log/i(o)) - /i(°)'^'=C/(i)^/i[°^ - ii7(i)°(log'/i(o)). (41) 



Together with Eqs. pip and ([40| . we can see that Eq. ([37| holds from the above. 

Furthermore, h^^jj satisfy the following equation [sec Eq. (I17p ] as a regularity condition 

- J(logM°))/^i;' - {h^'^^'h^^Ih^a^ + liDX'' + D.U^l (42) 
Then multiplying Ci to the above and using Eq. ()16p . we can see 

iA(i)A/i«=£/''"")i?,, + OCu-") (43) 
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holds. In the above, the higher-order terms like Cfh^^^, {Cih^^^Y , {Cih^^^Y , Cih^^^ Cih^^^ and so on are contained in 

0(u~"). Thus Eq. (|3T|) tells us that Cih^^^ = 0(m~") holds. As a consequence, we can show the m = 1 case of 
Eq. dSSl). 

For the m > 1 cases of Eq. (|33p , we perform the same procedure inductively. Let m > 1 be an integer and assume 
that the metric satisfies 



(Cr.fCigAB^o(\\ (44) 



for all fc(< to). Then {Cn)™CigAB =0{u ") is equivalent with 

=o(^\, (45) 



= O (L\ , (46) 



and 



'C'^aT = O - . (47) 



Since iJ^"*)" is written by fJ^^)", A^^), • • • , /i^™ from Eq. (P like Eq. the induction assumption 



^^/^S=0(^-j (z<TO-l) (48) 

immediately shows us that Eq. (|46|) holds. 

From Eq. ([7]), is written by U^^\ A^^\h!^2{i < m) and the trace part of h^'^^ ^ . Since the trace part of /i^j™"* 

is written by U^^\A^^^ and h'-^^ {i < to) through Eq. (gl), A^'"' is written only by U^^\ A'^'^^h'^^l {i < m) in the 

end. Then, by the assumption of the induction, Eq. p5|) holds. Next, the evolution equation for h''^^^ is described by 
Eq. Expanding Eq. ([6]) near the event horizon and multiplying, (£„)™~^, Eq. ([6|) becomes the following form on 
the event horizon 

h'^^ + \A'''^h^'=Fab, (49) 

where Fab is a function of h'^^l (i < to), h^^^ (i < to) and so on. Acting Ci to Eq. (|49|). then, we can see that 

(l/2)A(i)/i^™^ = Fab = 0(1/m") for j < m. Thus Eq. gT]) holds. 

Now we can confirm that the induction loop is closed. Then we can show that Eq. ([55]) . equivalently Eq. ([5^ holds 
if the spacctime is real analytic. Therefore we can show that the null geodesic generator / of the event horizon is the 
asymptotic Killing vector at late time in the sense of Eq. ([5^ . 



1 



IV. SUMMARY AND DISCUSSION 



We have confirmed that the expansion and shear of the event horizon should decay at late time in the vacuum 
spacetimes. Then, assuming the compactness of the cross sections of the event horizon, the null geodesic generators on 
the horizon give us an asymptotic Killing vector I at late time. This means that the horizon has late-time symmetry. 
By solving the Einstein equations, then, we have found that this late-time symmetry can be extended outside of the 
black holes. Therefore, at late time, there is the asymptotic symmetry outside of black holes. 

If the black hole rotates and there is another asymptotic Killing vector at late time, k, which will be a stationary 
Killing vector at a far distance and spacelike near the horizon, fc — / is also an asymptotic Killing vector expected to 
correspond to axisymmetry. In this sense, we would expect that the rigidity holds in gravitational collapse at late 
time. In these discussions, we assume the compactness of the event horizon. Thus this result cannot be applied to 
other null hypersurfaces which do not have a compact cross section. 

There is a remaining issue. In the proof of the rigidity theorem, the exact stationarity does show that the null 
geodesic generators of the horizon are Killing orbits. On the other hand, our argument could show us that the null 
geodesic generators of the horizon is a Killing orbit without assuming the presence of asymptotically stationary Killing 
vectors. It is likely that this difference suggests the existence of important and new points in black hole physics. 
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Appendix A: Einstein equations near event horizon 

In Appendix A, using the suitable variables, we will show the derivations of the Einstein equations near the event 
horizon. 

1. (n + 1) decomposition 

First we describe the formula of the {n + 1) decomposition. The metric can be written as 

gAB = enATlB +^AB, (Al) 

where jab is an n-dimensional induced metric, ua is the unit normal vector with uati^ — e — +1 {n^: spacclike) or 
— 1 (n^: timelike). 

We define the extrinsic curvature as 

Kab = \CnlAB- (A2) 

Now Ua can be written as ua = eN{d£L)A, where is a function which describes the hypersurface as = const, and 
N is the "lapse" function. Then the Riemann tensor is decomposed into 

RefghIa^Ib^Ic'^Id" = ^^^Rabcd - cKacKbd + cKadKbc, (A3) 
REFGDlA'^lB''lc^n^ = VaKbc - ^bKac, (A4) 

RACBon^n'' = -Cr^KAB + KacKb^ - e^VA^BN, (A5) 

where denotes the covariant derivative with respect to jab- 
The Ricci tensor is decomposed into 

RABU^n^ = -CnK - KabK'^'' - e^V^A^, (A6) 
Racti^Ib'' = ^^Kab - ^bK, (A7) 
RcdIa'^Ib" = ^^^Rab - ^CnKAB " ^KKab + 2eKAcKB'' - ^V^VsA^. (AS) 



The Ricci scalar is written as 



R - 2eCnK - eK^ - eKABK^^ - ^V^A^ 

+ eK^ - €KabK^^ - ^^^^ - 2e\/A{Kn^). 



(A9) 



The components of the Einstein tensor are 



GAsn^n^ = ^{-e^-'^R + - KabK""^), (AlO) 
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GAcn-^Ts'' = V^Xab - VbK, (All) 

GcdIA^Ib'' =^''^Gab - eKKAB + 2eKAcKB^ + ^-1ab{KcdK^'' + K^) 

1 1 (A12) 

- eCnKAB + eiABC,,K - —Va^bN + —-iab'^^N. 

2. Einstein equations 

We apply the (n + 1) decomposition presented in the previous seetion to the r-eonstant surface in our current four 
dimensional metric form: 

ds^ = -Adu^ + 2dudr + habidx" + U''du){dx'' + U^du). (A13) 
We express the above in the following form 

ds^ = N'^dr^ + g^^(da;^ + Nf'dr){dx'' + N^dr), (A14) 

where 

= J, (A15) 
= -j, (A16) 

= Lu". (A17) 
qfj,„ is the induced metric on the r-const. hypersurface as 

Note that the Latin indices a,b, ... and the Greek indices fJ-jV,.. are raised and lowered by hab and q^^, respectively. 
The unit normal vector to the r-const. hypersurface is given by ttia — N[dr)A and = N~^{dr — N^^d^)^. 
The extrinsic curvature on the r-const. hypersurface is defined as 

K^, - ]^C^q^, = -^{drq^j,,-V^N,-V,N^), (A19) 

where 2?^ is the covariant derivative with respect to q^^. 
We rewrite the induced metric as 

qf.udx'^dx" = -a^du^ + habidx" + (5''du){dx'' + P^du), (A20) 

where 

= A, = U". (A21) 

The timelike unit vector to the u-const. surface is given by m = —adu and = a~^{du — I3°'da)'^. 
The extrinsic curvature on the u-const. surface is given by 

kab - ]:C.uKb = ^{duhab-DaPb-DbPa), (A22) 

where Da is the covariant derivative with respect to hab. 
We introduce the following quantities for later convenience 

e = Ku,u^'u'' = -la,(loga) - log A^, (A23) 
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p"^ = = -drfi'' + log a, (A24) 



Kab = Kcdha^ht'^ = ^drhab + kab, (A25) 

and 

K = Kabh"' = ^{logh)' +k, (A26) 

where p^^u^ = = Kf^^u^, h = det hat and k = kath"'''- The prime denotes the r differentiation. Then Kf^^ can be 
written as 

K^^ = Qu^Ui, -2p(f,u^) + K^j,^. (A27) 

Using these quantities, we can decompose the four dimensional Ricci tensor Rab into the quantities on two dimen- 
sional space: 

RABm^m^ = -i(e (6 -K)-e2 + 2p,p"-K,bK"'' 

+ ^{CuCuN + kCuN - D^N - D^NDa log a), (A28) 

RABm^q^cu^ = -CuK + D'^Pa - kab^^'"'' + "^P^Da log a - Qk, (A29) 

RABQ^cq'^DU^u'' = -CuS + e^^ -eK-2p''pa-2p''Da\og — 

N a 

+ D^logaDalogN + -D^a ~ Cuk ~ kabk^'' ~ ^CuCuN, (A30) 
a N 

RABm'^h^'' = eL»"loga-2p6r''-fcp" + L)faK"''-L)°K + 5''e + K"''56loga--(a„p''-£«p"), (A31) 

a 

RABq^cu'^h^'' = Dak'''^D''k^p''K-2paK''''-^D'CuN + t^''Dal0gN-^ip')' 

- QD" log - - n'^^Da log - - -{dup" - /I/jp"), (A32) 
a a a 



and 



RABh'^ah^b = ^'^^Rab + Cukab + kkab - 2kackb'' - —DaDba - -^K^b - —kab + —CpHab - 2p(aDb) log — 

a N a a ^ ' a 



+ (e - K)Kab - 2paPb + 2HacHb'' - ^DaDbN + kabCa log N, (ASS) 



where the dot denotes 9„. The vacuum Einstein equation is given by Rab — 0. 

Appendix B: Explicit form of Einstein equations 

In Appendix B, we describe the components of the Einstein equations in terms of the metric functions explicitly. 
Using 

u = a-\l-U''da)=a-\du-U''da), (Bl) 
e = -(Ai/2)' + i£„logA (B2) 
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= ^f/"' + ^^"logA, (B3) 

and 

, 1 . _ 

l^ab = -^^ab + ^^JTJ^^^^b - DaUb - D^Ua), (B4) 

we can rewrite the Einstein equation Rab = in terms of our metric form. We will not provide all components of 
the Einstein equation explicitly. 

For example, RAsh^ah^b = becomes 

CiK, + ^A'h',, - + jh^KcKd + \h'%Kchbd + K^hac) ~ \KciDbU' + D^u,) + CuK, 

- ^h^{DaU^ + D'Ua) - ^hachbdU^' W"' ~ ^[(log/l) - 2DaU'']h^, 

- ^{log h)'{Ah'^, + hab - DaUb - DbUa) - ^^fe + \{hbcDaU'' + hacDbU'')- (B5) 

This determines the evolutions of hab- For RAB'm'^h^'^ = 0, we have 

CiU"' = ^D^A' - h'^'-'U''' (hbc - DaUb - DbUa) + h'^^D^Ah,^ + he - DbUc - D,Ub) - AD''{\ogh)' 

- i[/"'[(log/i) - 2DaU''] - i(log;i)'.D"A-^"[(log;i) - 2DbU'']+CuU''' . (B6) 

Appendix C: The gauge issue for A'^' [Eq. (|40[)] 

In Appendix C we will show the presence of the gauge where Eq. (j40p is satisfied. In our coordinate x"^ ~ {u, r, a;°), 
the metric can be written as 

ds^ = -Adu^ + 2dudr + hab{dx°- + U°-du){dx^ + U''du). (CI) 
The components of the metric are expanded near the event horizon (r = 0) as 

^ = r^(i) +0(r2), (C2) 



and 



Here A^^^ can be decomposed as 



hab - h^a^ + 0{r). (C4) 
= A^^^ + A'^^\u,x''), (C5) 



where A^^'^ is set to be a constant as shown in Ref. [T^ . 

When we consider the gauge transformation x^ x^ + the metric is transformed as 

gAB gAB + C(^gAB = gAB + 5gAB- (C6) 

To keep our gauge, the following conditions will be imposed: 

Sgur = , bgrr = , ig^a = , 
^guu = 0{r) , 5gua - 0{r) , Sgab = 0(1). (C7) 

From Sgrr = 29^^" = 0, we have ^" = f{u,x°'). Since 5gra and 5gur are given by 

Sgra = daC + UadrC + habdr^\ (C8) 
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5gur = {-A + U''Ua)dre + drC + + UadrC. (C9) 

5gra = and Sg^r = lead to 

r = -rduf + daf Wdr', ^ = ~dbf f h'^'dr'. (CIO) 

Then Sguu becomes 

where we used Eqs. (jC2p and (jCSp . This means that A'-^^ is transformed under the gauge transformation as 

^ + 94/^(1)) + dlf. (C12) 

Thus if we choose / as 

/ = ^ / duA^'\ (C13) 



A^^^ is transformed as 



Ad) ^ = A« - ^[9ui('' + (CM) 

Let assume that A'^' decays as u — cx3. For the moment, we write A'^' = 0(l/ii'"), where m is an integer. If 
TO > n, Eq. ()40p is aheady satisfied. Therefore, we suppose that to is smaUer than n. 
In the current gauge transformation, the transformed A^^'f satisfies 

i^i' = 4'' +0(l/u'"+'). (C15) 
Repeating this procedure, we can always choose the gauge satisfying 



(CM) 



If one wishes, one can continue the same procedure and then achieve an arbitrary faster decaying rate. But, the above 
is enough for our current purpose. 
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